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The dynamics of a single impurity interacting with a many particle background is one of the 
central problems of condensed matter physics. Recent progress in ultracold atom experiments makes 
it possible to control this dynamics by coupling an artificial gauge field specifically to the impurity. 

In this paper, we consider a narrow toroidal trap in which a Fermi gas is interacting with a single 
atom. We show that an external magnetic field coupled to the impurity is a versatile tool to probe 
the impurity dynamics. Using Bethe Ansatz (BA) we calculate the eigenstates and corresponding 
energies exactly as a function of the flux through the trap. Adiabatic change of flux connects the 
ground state to excited states due to flux quantization. For repulsive interactions, the impurity 
disturbs the Fermi sea by dragging the fermions whose momentum matches the flux. This drag 
transfers momentum from the impurity to the background and increases the effective mass. The 
effective mass saturates to the total mass of the system for infinitely repulsive interactions. For 
attractive interactions, the drag again increases the effective mass which quickly saturates to twice 
the mass of a single particle as a dimer of the impurity and one fermion is formed. For excited states 
with momentum comparable to number of particles, effective mass shows a resonant behavior. We 
argue that standard tools in cold atom experiments can be used to test these predictions. 

PACS numbers: 


I. INTRODUCTION 

Ultracold atom systems are effectively used as a test 
bed for condensed matter models. They are preferred 
because of the high degree of control in experiments such 
as tunable interactions, impurities and direct measure¬ 
ments by optical techniques. Certain theoretical models 
of condensed matter such as resonant interactions [l| or 
bosonic Mott transition 0 have been realized for the first 
time using cold atoms. Many models of one dimensional 
systems have been realized usiM two dimensional optical 
lattices to form narrow tubes [j-Q. 

One of the powerful theoretical tools to describe one 
dimensional systems is the Bethe Ansatz (BA). BA so¬ 
lution has been generalized to many integrable models, 
e.g. systems with multiple components, different statis¬ 
tics or spin This exact solution method has been 

employed to explain experimental data on a number of 
instances [1,[HI. However, as BA methods are restricted 
to one dimension, they have not been used to describe 
systems where an external artificial gauge field is present. 

In one dimension, such an external magnetic field can 
be disregarded by using a gauge transformation, unless 
the one dimensional system closes onto itself. Thus, if 
the particles are confined to a ring as opposed to a tube, 
the artificial magnetic field will significantly effect the 
physics. Such rings, in the form of toroidal traps, have 
been realized experimentally [T^ - [l9| . Although none of 
these experiments have included an artificial gauge field 
so far. 

In this work, we consider such a toroidal trap contain¬ 
ing non-interacting fermions and describe the behavior 
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of a single charged impurity interacting with background 
atoms. We argue that an artificial magnetic field coupling 
to the impurity is an efficient way to probe the polaron 
state forming due to the interactions. Artificial magnetic 
fields are created W coupling light to the internal states 
of the atoms [20l - l^ . Hence, they are highly specific to 
the internal state making it possible to create effective 
magnetic fields coupling only to one type of atom. 

The charged particle is expected to drag the uncharged 
fermions along with itself around the ring. Because of the 
interactions between the impurity and the background 
atoms a collective excitation usually called a polaron is 
formed [l^. This excitation will couple to the external 
magnetic field with the charge of the impurity particle, 
however, its mass will critically depend on the interaction 
strength. The amount of angular momentum carried by 
the impurity and the uncharged fermions also depend on 
the total external flux through the ring. By changing the 
artificial magnetic field strength, it is possible to access 
excited states of the system adiabatically. We show that 
an artificial magnetic field coupling specifically to the 
impurity would be a very effective tool to probe polaron 
physics. 

We describe this system exactly using a Bethe Ansatz 
(BA) solution for contact interactions which are justi¬ 
fied for cold atoms as the dominant scattering is s-wave. 
For strongly attractive interactions, the impurity forms a 
bound state with one of the background fermions and the 
physics reduces to the motion of a dimer with twice the 
mass of the particle. In the other limit of infinitely re¬ 
pulsive interaction, effective mass saturates to total par¬ 
ticle number. We calculate the energy and momentum 
distributions, total transferred momentum and the effec¬ 
tive mass for all interaction strengths. We believe these 
results can be experimentally checked with state of the 
art toroidal traps and techniques for artificial gauge field 
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FIG. 1. (Color online) A simple illustration of the system. 
N — 1 uncharged fermions (light gray) and a single charged 
impurity (dark gray) are trapped on a ring. The impurity is 
interacting with the fermions via Delta-function interaction. 
An artificial magnetic field couples exclusively to the impu¬ 
rity. The dynamics of the system depends on the interaction 
strength between particles and the total flux through the ring 

P=^- 


generation. 

The paper is organized as follows. In the next section, 
we define the model, introduce the notation and review 
earlier studies. In Section III, we solve the system for 
two particles and then generalize to any particle number 
using the BA. Sec. IV contains the analytical solution of 
the BA equations in certain limits and comparison with 
numerical solutions. We present our results for several 
quantities such as energy, angular momentum and effec¬ 
tive mass of the charged particle. We give our conclusions 
along with a brief discussion of possible experiments in 
Sec. V. 


II. THE MODEL 


The first quantized Hamiltonian for one charged par¬ 
ticle among N — 1 uncharged fermions under a magnetic 
field reads 


\ (h d n2 fi? ^ ^ 

j=2 3 j=2 


( 1 ) 

All particles are assumed to be on a ring of radius R, 
0 < Xi < 2 'kR. The position of the charged particle is 
xi and A is the vector potential in the symmetric gauge. 
The Hamiltonian can be made dimensionless by using, 
E = d = and /3 = /3 is 

the total magnetic flux through the ring in units of flux 
quantum q/h. Dropping the tildes 


f) 2 ^ rP" ^ 

j=2 3 j=2 

The effect of the magnetic field can be shifted to the 
boundary conditions by a gauge transformation [^ . 
Namely, when the first particle makes a full circle around 
the ring the wave function gains a phase factor of 


where the periodic boundary conditions (PBCs) for the 
uncharged particles remain unaffected by the gauging 
process, 

n ^ (3) 

Apart from the twisted BCs, the 5-function interaction 
can be handled as a two-sided boundary condition (BC) 
between two different regions of N-particle space corre¬ 
sponding to different permutations of particles. The dis¬ 
continuity relation at the boundary Xi = Xj (which is 
obtained by passing to the center of mass and relative 
coordinates and then integrating the Hamiltonian) is 




X\ <Xj 


(dj-di)^ 


XjKxi 



j ^ 1- 


( 4 ) 

This one dimensional problem of two-component 
fermions has been studied by using the BA in the pre¬ 
vious century. First, the one-spin deviate problem in a 
Fermi sea is solved by McGuire @ and Flicker and Lieb 
solved the two-spin deviate problem. Yang Q ele¬ 
gantly derived the BA equations for the general M down- 
spins among N up-spins. Twisted BCs have been used 
throughout the BA literature as a way to probe ground 
state properties. However, with the possibility of opti¬ 
cally inducing artificial magnetic fields, it is important 
to calculate the properties of the system at finite flux as 
opposed to infinitesimal values near zero. It is also neces¬ 
sary to consider cases where different components in the 
system experience different gauge fields. 

Our calculation takes both of these constraints into 
account and allows us to exactly study the dynamics re¬ 
sulting from the dragging effect of the charged particle 
on the uncharged particles. The resulting polaron physics 
has attracted great interest in the context of cold atoms 
over the last few years [ 21 , 123 . 


III. THE ANSATZ 

As the interactions are reduced to BCs, the wave func¬ 
tion for a given permutation of the particles is a super¬ 
position of plane waves. As the collisions of equal mass 
particles in one dimension conserve magnitudes of the in¬ 
coming momenta, the interacting problem is integrable. 
Hence, in a given region only a finite number of plane 
waves are needed to construct the wave function. To 
make our notation clear, we first start with the case of 
one charged particle with one neutral particle. 


A. N=2 Particles 


For two particles we have 2! = 2 regions and the wave 
function in these regions is expressed as follows: 


'^i2{xi,X2) = -k 

4'2i(a;i,X2) = (12)2ie*('=i“i+'==^^) -k (21)2ie*('==^^i+'=^^^)(5) 
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where we use parenthesis with a subscript to indicate the 
coefficients of plane waves. In this notation numbers in 
the parenthesis indicate the order the wave vectors fci, ^2 
are distributed to the coordinates in the exponent and the 
subscript indices indicate the ordering of the particles on 
the ring, i.e. 'hi2 means xi < X2- At xi = X2, the wave 
functions in the two regions should be equal whereas their 
derivative should obey Eq 01 Equating the coefficients of 
each plane wave on both sides, we obtain: 

BCs: at Xi = X2, 


( 6 ) 


Sl 2 A 

(7) 


( 8 ) 

Allowed values for ki, k2 are found by applying the PBCs. 
PBC for one of the particles gives the BA equation. 
BCs: at 27 r 

as a:2 : 0 ^ 27 r, 'I'2i(a;2 = 0) = 'Si2{x2 = 27 r), 

(12)21 = (12)126*'=="^, (21)21 = (21)i2e*'=i2-. 

(9) 

as a;i : 0 ^ 27 r, 'I'i2(a;i = 0 ) = = 27 r), 

(12)i 2 = (12)2ie*('=i+^)2-, (21 )i 2 = (21)2ie*('==+'^)2-. 

( 10 ) 

Combining the two BCs at 2 tt, we obtain another con¬ 
straint fci -b ^2 -|- /3 = n, for n S Z. This is a reflection of 
the total angular momentum conservation in the system. 
Eqsl 5 ]and EqlTUlhave non-trivial solutions only when the 
determinant below vanishes, 


(12)21 -b (21)21 = (12)i2-b(21)i2, 

(12)21 — (21)21 = (12)12^1 + -—^ + (21)12^ — 1 

where S12 = *(fci — k2)/c. Combined BCs give 

/(12)\ _ + 7 ^ \ Al2)\ 

U)j2i"Vili i-^/Wi 


j_1_g- 2 (fei+/ 3 ) 27 r 

~ S 12 


1 

S 12 


= 0 . 


-1 \ _ 1 _g- 2 (fc 2 +^) 27 r 

S12 S12 

( 11 ) 

Solution of this determinant gives the BA equation, 
a = |cot (^|(a-bn-/ 3 )) (|(«-'«+ /?))> (12) 

where energy is A = fQj- a = k2 — ki. For the 

two-particle case, this problem can also be solved exactly 
without using the BA [? 7 l |. 


'cos(7r(n — /?)) 
. sin(7rQ:) 


— cot( 7 rQ;)^. 


(13) 


These two equations analytically reproduce each other 
and the numerical results match perfectly (Figl 2 |). 

Extension of this method to N particles is straightfor¬ 
ward if cumbersome. 



FIG. 2 . (Color online) Energy of the lowest three states 
vs. interaction strength for N — 2 particles, for zero to¬ 
tal angular momentum. Only scattering states are displayed. 
Energy is calculated by three different methods. Lines are 
from Eg 1 131 direct analytical solution without employing BA, 
which is algebraically same with the two-particle BA calcula¬ 
tion lEa ll2l) . Diamonds are from the general A'-particle BA 
calculation ('Ea l 21 ll . 
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FIG. 3 . (Color online) Ground state energy vs. flux P for 
N = 2 particles with total angular momentum n. Eigenstates 
for flux (3 with total angular momentum n are also eigenstates 
for flux / 9 -bl with total angular momentum n-bl. The system 
can be analyzed by considering flux values between —1/2 < 
/3 < 1/2 for all n. As the flux is increased by one adiabatically, 
the system evolves to a higher excited state which has one 
more unit of total angular momentum. As can be observed, 
the crossings between different eigenstates is not a problem 
for adiabatic evolution since only states with different total 
angular momentum n are degenerate. 
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B. N-1 Fermions, One Charged Particle 


The distinguishable charged particle is denoted again 
by xi and the wave function is defined in TV! regions 

I 


corresponding to different permutations In each one 
of these regions the wave function consists of Nl plane 
waves in its most general form without imposing the an¬ 
tisymmetry between the fermions. As a total we have 
N\ X N\ coefficients: 


123.. 

. = (123. 

\ i{kiXi+k2X2+k3X3. + ... 

• •) 123... 

) -b (213. 

• ■)l23.. 

^i(k2X-i_+k3X2 + k3X3 + .. 

..) + 

213.. 

. = (123. 

\ i{klXl+k2X2 + k3X3 + ... 

■ •j213...e 

) -b (213. 

■ •)213.. 

^i(k2Xl+klX2 + k3X3 + .. 

..) + 

132.. 

. = (123. 

\ i{k3X3+k2X2 + k3X3 + ... 

• •) 132... 

) -b (213. 

• ■)l32.. 

^i(k2X3+k3X2 + k3X3 + .. 

..) + 


(14) 


r 


where ki, k 2 , ■ ■ ■ k]\f are distinct wavenumbers. BCs at 
Xi = X 2 are not effected by the addition of other fermions 
at the end of the sequence: 




213... 

BCs at 27r follow the same logic; 

as a :2 : 0 27r, 4'2i3...jv(a:2 = 0) = 'i>i 3 ...N 2 {x 2 = 2tt), 

yielding 

(123 . . .)213...W = (123 .. .)i3...Ar2e**^^’^, 

(213 .. .)2i3...Ar = (213 .. .)i3...JV2e*''^"". (16) 

Number of independent coefficients decreases consid¬ 
erably by requiring antisymmetry upon exchange of 
fermions. Every coefficient of a plane wave in region 
xi < xs < X 2 < ... < xn is identical with the co- 
efhcient of the same plane wave in region xi < X 2 < 
X 3 < ... < Xn- This can be shown by noticing that at 
X 2 = X 3 the wave functions must vanish requiring e.g. 
(123 .. .)i 23 ...Ar = —(132 .. .)i 23 ...Ar. Fermionic antisym¬ 
metry also relates the wave function in separate regions 
4^123...IV = —4^132...TV. As a result, the coefficients only 
depend on the position of the charged particle in the or¬ 
der. We can move indistinguishable fermions through one 
another at will and the N\ regions reduce to N regions. 

After this simplification it is easy to combine the BC 
at a 5-function with the overall PBC. 



Sl2 Sl2 


^(123...)^ 


213...iV 


-1 

Sl2 


l-jfe/ V(213...), 


satisfy 

c c c 

/ci — - cot Tvki = k 2 — — cot TTk 2 = ks-- cot irk^, = ... = A, 

where A is a real constant. This form is equivalent to the 
usual BA equations Q- 

Hence, the N wavenumbers which define an eigenstate 
must be chosen as N distinct roots of the equation: 


k — X = - cot nk. 
2 


(19) 


However, there is another constraint. Applying PBCs 
sequentially on all particles restricts A. 

As xi : 0 ^ 27r, 4'i23...Ar(a:i = 0) = e*^^’"4'23...Ari(a:i = 27r), 
(123)123 = (123)23le*(^+'=^)'^ 

as a ;2 : 0 27r, ( 123)231 = ( 128 ) 3126 *^^^^, 
as X 3 : 0 ^ 2 ^, (123)312 = ( 123 )i 23 e*'=^ 2 -, 

In combination: 

(123...)i 23 Tv = e*('=i+'==+-+'=^+'^)""(123...Af)i23 tv 

( 20 ) 

reflecting angular momentum conservation, sum of all the 
wavenumbers plus the flux must be integer on a ring. 

In short, the BA equation is solved by finding N roots 
of a simple equation subject to the angular momentum 
constraint: 


123...TV 


V'==2-(123...)i23...TV^ 

^e*'=i2-(213...)i23...TVy 
The determinant can only vanish if ki and k 2 satisfy, 
c c 

ki — -cotTrfci = k 2 — -cot7rfc2. (18) 

The same procedure can be applied to any pair of 
wavenumbers ki,kj. Thus, all the wavenumbers must 


(17) 


k — X — — cot irk, 
2 


TV 


kj = n — j3, 


n e Z. (21) 


So far our treatment implicitly assumed repulsive interac¬ 
tions. In which case, all the wavevectors ki are real. The 
ansatz can easily be extended to attractive interactions 
yielding exactly the same equations EqUHjli]. However, 
for negative c two of the roots will be complex, as the 5- 
potential in one dimension has only a single bound state. 
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FIG. 4. Ground state energy vs. interaction strength for 
N = 1000 particles and zero total angular momentum. Nu¬ 
merical solution of BA equation (dots) Ea l21l are virtually in¬ 
distinguishable from analytical solution (circles) E = (Fermi 
energy of — 1 fermions)-|-Ai?. Error between the numerical 
and analytical solutions are too small to observe even in the 
regimes where the assumptions for analytical calculation fails. 

IV. SOLUTION OF THE BA EQUATION 

The cot nk term in Ea l21l diverges at every integer k, 
thus, regardless of the value of /3 (or A) there is a root 
between every consecutive integer (Figl^. By changing 
the value of A, all roots can be adjusted so that the to¬ 
tal angular momentum constraint is satisfied. All the 
eigenstates in this problem can be labeled identically by 
choosing N distinct integers corresponding to the differ¬ 
ent branches of cot nk and the total angular momentum 
n G Z. The energy of an eigenstate is simply the sum of 
squares of all wavenumbers 

N 

E = Y^kl ( 22 ) 

i=l 

For the simplest case of /3 = 0, the ground state corre¬ 
sponds to A = 0 and the total angular momentum n = 0. 
The roots k are distributed symmetrically around zero 
for even N, hence, automatically satisfy the total angular 
momentum condition. The wavevectors for the ground 
state are in the N branches of cot from —N/2 to N/2. 
Excitations above this ground state can be generated by 
two procedures. First, by changing A, N roots which are 
on the same branches of cot can be generated so as to 
create an eigenfunction with non-zero total angular mo¬ 
mentum (n ^ 0). Second, at least one of the roots can 
be chosen to reside on a branch that is not occupied for 
the ground state. For such a particle-hole excitation A 
must be adjusted to ensure the total angular momentum 
constraint. 

Inclusion of the magnetic field affects only the total 
angular momentum constraint. As that constraint is de¬ 



FIG. 5. (Golor online) A representation of graphical solution 
to BA equation. The cot nk term in Ea l21l diverges at every 
integer k and there is a root between every consecutive integer 
independently from A(/I). By changing the value of A, all 
roots can be adjusted so that the total angular momentum 
constraint is satisfied. 


fined only up to an integer (n), the problems with values 
of j3 differing by an integer are identical. Eigenstates for 
flux j3 which have total angular momentum n are also 
eigenstates for flux ,8-1-1 which have total angular mo¬ 
mentum n-|-l. This is a restatement of flux quantization. 
We can analyze the system by considering flux values be¬ 
tween — 1/2 < P < 1/2. 

However, in an experimental setting slowly increasing 
the value of the flux through the ring is a useful method 
to access excited states. As the flux is increased adia- 
batically from zero to one, the ground state evolves to 
an eigenstate which has its roots exactly in the same 
branches as the ground state, but, has a total angular 
momentum of minus one at zero flux (FiglH]). The cross¬ 
ings between different eigenstates do not pose a problem 
for adiabatic evolution as only states with different total 
angular momentum n can be degenerate in energy. 

The BA equation Eql^can be very efficiently solved 
once the regions for the roots are determined. We used 
the Newton-Raphson algorithm to find a solution within 
a particular region. As all the roots depend monotoni- 
cally on A, another Newton-Raphson search is employed 
to satisfy the total angular momentum condition. We 
have found numerical solutions for systems of up to 10000 
particles with high accuracy. 

Although numerically solving the BA equation is effi¬ 
cient and accurate, an analytic solution can provide more 
insight about the physics of the system. Analytic formu¬ 
lae for energy, angular momentum and effective mass also 
would be desirable to make correspondence with experi¬ 
mental observations. 

In the limit of strong interactions 1/c -C 1 and large 
particle number N/c^ 1, such an analytic form can be 
obtained by approximating the roots of the BA equa- 
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tion. In this limit, because the cot diverges quickly near 
integers, most of the roots are close to integers. Apart 
from the few roots near k ~ Ajthe deviation of the root 
A from an integer s is small Q- Solving for this small 
deviation we find that the roots occur at 


is small compared to N and the integral can be approxi¬ 
mated as 


n — P 


c 


pXF+b 

J xp — b 


dx atanx 


cb 

—atanxF. 

TT 


(25) 


kt=s H—acot-(s — A), 

TT c 

12 N 

k7=—s -acot-(s -I- A), s = 0,1,...,--—1, (23) 

TT c 2 


Through this relation b, hence A, is obtained for any flux 
value, allowing us to find expressions for all the roots in 
a self consistent way. 


where acot is defined in the continuous region (0,7r) for 
EqsHnito be accurate guesses. Here we have restricted s 
to analyze the ground state and excited states with roots 
on the same cot branches. Applying the total angular 
momentum condition we get. 


N/2-1 

n - /3 = ^ ks 

s=0 


1 

TT 



C 

27r 



da; ^acot(a; — b) — acot(a: + &)^, (24) 


with b = 2\lc and xp = {N — Xj^^jc. Here the initial 
assumption of strong interactions and large particle num¬ 
bers allow us to approximate the sum by an integral. For 
the ground state and the first few excited states n — p 

I 


A. Energy 

Using these expressions for the roots, the total energy 
is 

Ar/2-l 

E=Y^ kl 

s^O 

f 2 2s / 2 2 \ 

= s 2s -I -( acot-(s — A) -I- acot-(s -I- A) 1 

^ TT C C 

-l-^((acot-(s - X)Y + (acot-(s -|- A))^^|. (26) 

The first term above is the total ground state energy of 
N — \ non-interacting fermions. Interactions result in the 
second and third terms which are first and second order 
corrections in our expansion. When A’s are small, the 
third term is negligible. In this limit, the energy shift 
due to interactions is 


for repulsive interactions: 
AA = ch{n -P) + ^ 


^(x_F -I- 6)^ -I- l^atan(a;F + 6) + {{xp — b^ + l^atan(a:F — b) — 2xp 


(27) 


with 


7r(n — P) 
c atan(xF) 


This approximate form for energy successfully repro¬ 
duces numeric results for particle numbers as small as 4 
throughout all the interaction range. Ground state en¬ 
ergy as a function of interaction strength is plotted for 
a typical case in FiglUfor 1000 particles at P = 0.2 flux. 
The deviation between numerical and analytical results 
are too small to observe in this plot. 

For attractive interactions, the (5-function interaction 
supports one bound state in one dimension. Correspond- 

I 


ing imaginary wavevectors appear as solutions of the BA 
equation. For k = a + ia with (a, a) G K, the BA equa¬ 
tion has only two roots with a ^ 0. The charged parti¬ 
cle is bound with only one of the background fermions. 
When l/|c| <C 1, the complex roots are at A: = A ± zc/2 
while the rest of the roots preserve their form of Fasl23l 
If the bound state is narrow, Pauli repulsion between 
the fermion in the bound pair and background fermions 
decreases the effective interaction. Within these approx¬ 
imations, we analytically calculate the total energy. 


for attractive interactions: 

AE = ^ ^ —--|- cb{n + ^(a;_F-|-&)^-l-l^atan(a;F + 6) -I- ^(a:_F — 6)^-1-atan(a:_F — b) — 2xp\, (28) 
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with b = 


in- P) 


c^l — iatan(a;F)^ 



FIG. 6. Angular momentum of the impurity vs. interac¬ 
tion strength for N = 100 particles for (a) attractive, (b) 
repulsive interactions from the analytic calculation, numeri¬ 
cal solutions produce the same results. In the non-interacting 
limit, the impurity carries all the angular momentum (n —/3). 
Li saturates to almost zero for infinitely strong repulsive in¬ 
teractions as the total angular momentum is shared equally 
between all particles. The same behavior holds for excited 
states. For strongly attractive interactions, Li saturates to 
half the total value signifying dimer formation with one back¬ 
ground particle. The insets in both figures focus on the weak 
interaction limit. 


this system we are interested in how the dynamics of the 
impurity particle is affected by the fermion background. 
To this end, it is instructive to calculate angular momen¬ 
tum carried by the impurity Li and the related effective 
mass. As the impurity is interacting with the fermions, 
this effective mass is not only the mass of the impurity 
but also gets a contribution from the fermions dragged 
along with it. Such a compound object is generally called 
the polaron state or dimer state especially for attractive 
interactions. 

As stated above, one of the most interesting physical 
quantities in this system is the angular momentum car¬ 
ried by the charged particle, represented by the operator 
Pi = —As this particle is coupled to the exter¬ 
nal magnetic field, Pi is the canonical momentum not 
the kinetic momentum. However, canonical momentum 
is the quantity that is generally measured by expansion 
imaging in artificial magnetic field experiments. The ex¬ 
pectation value of {Pi) = Li is easily obtained by taking 
the derivative of the total energy with respect to flux. 


B. Angular Momentum 

To understand the physics of the system and make cor¬ 
respondence to possible experiments, it is important to 
calculate other measurable quantities. In particular, for 

I 


_ -1 d/\E 
5/3 


(29) 


Using the approximate form for the energy Ea ]27l we 
obtain 


7r(n — /3) 
atanx p 


iatanxp 


{xp + b)ataii{xp + b) — {xp — 6)atan(a;_F — b)>. 


(30) 


This form is valid for positive c and easy to interpret. 
In the non-interacting limit, the canonical momentum of 
the charged particle is fixed by the external flux. Hence, 
all the angular momentum is carried by the charged par¬ 
ticle. As interaction is turned on, the charged particle 
drags the background fermions and transfers some of its 
angular momentum to them. Stronger interactions in¬ 
crease the fraction of the transferred angular momentum 
and in the limit of infinitely repulsive interactions angu¬ 
lar momentum is equally shared by N particles. On the 
other hand, for strong attractive interactions, Li satu¬ 
rates to half of the angular momentum in the system 
proving the formation of a dimer with one background 
fermion. 

The behavior of Li is displayed in FiglS] and FiglT] 
as a function of c and /3. Even for N = 100 parti¬ 


cles, the difference between numerical calculation of the 
derivative and the expression given above is negligible. 
As a function of interaction strength, the rapid decrease 
and eventual saturation of Li validates the scenario dis¬ 
cussed above. The linear dependence on flux is expected, 
however, the slope of Li decreases as interaction gets 
stronger. This slope carries valuable information as it is 
related to the effective mass of the composite excitation 
formed by the impurity and background fermions. 

C. Effective Mass 
We define the effective mass as 




(31) 
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FIG. 7. (Color online) The angular momentum of the charged 
particle vs. flux at varying interaction strength for N — 100 
particles. As expected Li depends linearly on flux. The slope 
of the line decreases with increasing interaction strength in¬ 
dicating higher values of effective mass of the impurity. 
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-p=0.2 

N=50 


--- p=5.2 

n=0 


. p=10.2 


FIG. 8. (Color online) Effective mass of the impurity vs. in¬ 
teraction strength for = 50 particles and zero total angular 
momentum n = 0. (a) For attractive interactions, m* satu¬ 
rates to twice the mass of the impurity due to the formation 
of a tightly bound pair. The inset shows the behavior around 
zero interaction in more detail. For attractive interaction, the 
effective mass (given by Ea l32ll is almost insensitive to flux 
change, (b) Eor repulsive interactions, m* converges to N. As 
flux P increases, this saturation gets faster. The dependence 
on the flux is more prominent for small particle numbers. 


In the non-interacting limit, the effective mass is equal 
to m, however, its behavior is very different for attrac¬ 
tive and repulsive interactions. As repulsive interactions 
are increased, it gets harder for the impurity to tunnel 
through the fermions and the dragged particles increase 
the effective mass. The increase in the effective mass 
saturates only when all the particles are moving together 
with the impurity. Thus, at large repulsive interaction, 
the effective mass reaches Nm. For weak attractive inter¬ 
actions, the first effect is once again the drag increasing 
the effective mass. However, the attractive (5-function 
has a single bound state in one dimension. Thus, the im¬ 
purity captures one of the background fermions and as 
the size of the bound state gets smaller, Pauli exclusion 
effectively repels the other fermions. The effective mass 
for attractive interactions increases and reaches 2m for 
infinitely attractive interaction where a dimer is formed 
from the impurity and one fermion. For attractive inter¬ 
actions, the analytical expression in the strongly inter¬ 
acting limit is useful to calculate the dimer mass. 


5/32 
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We calculated the effective mass numerically and analyt¬ 
ically. For the ground state, the above scenario is vali¬ 
dated by these calculations (see FiglS]). The dependence 
of the effective mass on the external magnetic field is 
strongest for small particle number as this limit is the 
strongly interacting limit in one dimension. As the num¬ 
ber of fermions increases, effective mass in the ground 


state has weak dependence on /3. In this case, effective 
mass is essentially determined locally as the ability for 
the impurity particle to complete a full rotation is ham¬ 
pered. 

The utility of an external magnetic field is the access 
it provides to excited states through adiabatic pumping. 
Excited states in this system are expected to be stable 
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FIG. 9. Resonant behavior in m* for j3 comparable to 
N. When the drag effect applied by the background parti¬ 
cles overcomes the driving force of the magnetic field, m* can 
become negative. When the second derivative of the energy 
with respect to flux becomes zero, m* diverges. This diver¬ 
gence does not change the infinitely strong interaction limit. 


due to angular momentum conservation. It is thus rea¬ 
sonable to expect effective mass measurements to be car¬ 
ried out on such states in a cold atom setting. For the 
excited states, with angular momentum |n| < iV/4 the 
main effect is faster saturation of effective mass as c in¬ 
creases. However, for higher excited states, there is reso¬ 
nant behavior (Figl^. Due to the nature of BA solution, 
a state for which all the roots are on the cot branches 
from —N/2 to N/2 can have at most N/2 units of an¬ 
gular momentum. When the total angular momentum 
of an excited state is comparable to particle number N, 
the sharing of this angular momentum between the im¬ 
purity and the background is limited by this constraint. 
Thus, it is possible for this system to support negative 
effective mass if the external force acted by the magnetic 
field is overcome by the back reaction from the fermions. 
We numerically find this behavior for both low and high 
particle number (see FiglHl). Experimentally this effect 
should be more accessible for small number of particles 
as it is easier to pump angular momentum comparable 
to particle number. 


D. Correlations 

Apart from the single particle properties related to the 
impurity, it is instructive to look at global properties to 
understand how the external particle disturbs the one di¬ 
mensional Fermi liquid. A common way to visualize the 
disturbance in the Fermi sea is to plot the deviation of the 
distance between the BA roots (wavevectors) from one. 
For an undisturbed Fermi sea, this deviation is always 
one. For a weakly interacting impurity, the deviation 
is confined to a narrow region in k-space around n — P 
(see Fig ITU)) . This is expected as the impurity carrying 


1.08 


1.06 


1.04 

I 


1.02 


1 


"■'-50 -40 -30 -20 -10 0 10 20 30 40 50 

( k .+ kj^,)/2 

FIG. 10. Effective momentum density in k-space for different 
interaction strengths. In the strongly interacting limit, the 
distance between adjacent BA roots (wavevectors) is one. For 
small c, the roots are closer to each other around (n — fi). 
The impurity carrying n — /3 units of angular momentum in 
the non-interacting limit first disturbs the fermions which are 
momentum matched to that value. 



n — P units of angular momentum in the non-interacting 
limit first starts dragging fermions which are matched in 
momentum. As the strength of repulsion increases, so 
does the effected region in k-space, however, the devi¬ 
ation gets smaller. For infinitely repulsive interactions, 
the impurity becomes indistinguishable from the back¬ 
ground fermions iFig fTO)) . For highly excited states where 
n is comparable to N, particle-hole excitations complicate 
this picture similar to the effect we discussed for the ef¬ 
fective mass. 

Another important physical property is the two- 
particle correlation function. Although for d-function in¬ 
teractions only the value of this function at zero deter¬ 
mines the interaction energy, its general form is exper¬ 
imentally accessible through Hanburv-Brown-Twiss[^ 
type measurements. This correlation also can be re¬ 
garded as the real-space form of the bound state cre¬ 
ated by the impurity. To calculate the two-particle cor¬ 
relation function, we need to determine the coefficients 
of the plane waves in each region. Following McGuire 
we choose the first coefficient in the first region 
xi < a;2 < ... < xn , 

(123..iV)i23..iv = (1 - (33) 

Other coefficients in this region determined by BCs yield 
very similar expressions. The wavenumber associated 
with the distinguishable particle appears in the exponent 
and the sign of the permutation multiplies the coefficient: 

(213..iV)i23..iv = -(l-e*'"'==), (34) 

(312..Af)i23..jv = (35) 
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FIG. 11. (Color online) Two-particle correlation function for 
N = 2 particles. As expected, 312 at zero separation decreases 
with increasing interaction strength. For weak interactions, 
the correlation function at zero decreases with increasing flux. 
However, for strong interactions, the correlation is almost 
insensitive to flux change due to the fermionization of the 
charged particle. 


The coefficients in other regions are related to the same 
coefficient in the first region with a phase factor deter¬ 
mined by the PBCs. This phase factor is a full circle 
rotation around the ring of the particles that Xi has to 
pass to be in the given region. Thus, the momenta be¬ 
longing to the particles that xi has passed multiply the 
coefficient, e.g. 

(21354..Af)23i4..iv = (1 - ... (36) 

In the simple form, the wave functions are not normal¬ 
ized, but we normalize the correlation function at the 
end. The two-particle correlation function in any state is 
given as 


P‘2'K p27T 

gi2{xi,X2)= / dxs-- dxN'd'*'^'- (37) 

Jo Jo 

For the N = 2 particle case, the correlation function is 
simply the absolute square of the wave function. As could 
be expected, the correlation function is highly affected by 
the flux for the two-particle case. Using the numerically 
and analytically found wavenumbers in the expression, 

gi 2 ixi,X 2 ) = 4 — 2cos(27rfci) — 2cos(27rfc2) 

— 8 sin( 7 rA:i) sin( 7 rfc 2 ) cos(/ci — k 2 )(x 2 — Xi — tt), 

(38) 

we observe that inclusion of the flux generally decreases 
the two-particle correlation function (Fig fTTl) . However, 
if the interactions are strong enough so that the particles 
are almost fermionized, this decrease is very small. It is 
also notable that although the flux breaks time-reversal 
symmetry and the wave functions choose a direction on 
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FIG. 12. (Golor online) Two-particle correlation function for 
A’ = 50 particles, (a) For weak interactions, Friedel oscil¬ 
lations occur as interference of two waves with wavelengths 
related to kp — d and kp -\- jS. (b) At strong interactions, the 
correlation becomes zero at zero separation since the impurity 
is effectively indistinguishable and gi 2 and the frequency of 
Friedel oscillations are almost insensitive to flux change. 


the ring, correlation function is even with respect to xi — 
X 2 - This property holds for any particle number. 

For the general N particle case, we arrange the wave 
function in a better form to evaluate the integrals. We 
assign two wavenumbers to xi and X 2 , and the rest of the 
particles are represented by a Slater determinant since 
they are indistinguishable fermions. For example, if we 
have fci, /c 4 associated with xi^X 2 respectively, the Slater 
determinant is represented by D 14 indicating the use of 
all wavenumbers except ki and k^ in the exponents, 




^ik^Xi ^iksxi ^ik^Xi 


klA = 


(39) 
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Hence, the wave function in the first region can be written 
as, 

'h =^(12 ... iV) 12 ... AT + (21 . . . iV)i2...iV^Dl2 

+ ^(13 . . . iV)i2...1v + (31 . . . iV)l2...7V^Di3 + ... 

(40) 

Integrating 'I'*'h over X3 ,..., xn, the Slater determinants 
are orthogonal in large particle number limit, as the outer 
roots of cot’s are very close to integer values. The corre¬ 
lation function is then expressed as a sum over pairs of 
momenta associated with Xi and X 2 , 

N N 

(^) = E EI I' +1 I' (41) 

t<S 

-f 2i?e{(ts ... N)*{st... (42) 

where x = X 2 — Xi. Here, x > 0 but the symmetry of 
the correlation function for x < 0 can be easily seen by 
using the region X2 < X3 ... xn < xi instead of the first 
region. Finally, the correlation function is normalized to 
average density on the ring so that it saturates to one. 
The correlation function for iV = 2 and = 50 particles 
are given in FiglTT] and FigH^ respectively, for different 
interaction strengths and flux. As expected, the corre¬ 
lation function at zero separation 512(0) decreases with 
increasing interaction strength until it saturates to zero. 
The other important feature of the correlation function 
is the Friedel oscillations [1^ reflecting the sharpness of 
the Fermi surface in one dimension. 

The two-particle correlation function is a local quan¬ 
tity while the external flux changes the system properties 
globally. For any pair to feel the effect of the flux, one of 
the particles must go a full circle through the ring. Hence, 
as could be expected, the effect of the artificial magnetic 
field on the correlation function decreases as the number 
of particles increases or if they interact strongly. Even 
for the lowest lying excited states, we find the primary 
effect of the flux is on the Friedel oscillations while the 
shape of the correlation hole is unchanged. 

Finally, we calculate a thermodynamic quantity which 
is also related to 512(0). Derivative of the energy with 
respect to interaction strength c gives us interaction po¬ 
tential, so, the kinetic and interaction contributions to 
the total energy can be separated. In FiglT^l one can 
notice that the interaction energy makes a peak and 
then decreases for increasing c. The initial increase is 
expected, however, as interactions become stronger, the 
tendency of the fermions to avoid the impurity dominates 
and the impurity is effectively fermionized. This is appar¬ 
ent in the d-function BC EqHl Additionally, the interac¬ 
tion potential is equal to the correlation function at zero 
512 ((2^2 —2:1) =0) times the interaction strength which 
reproduces the results obtained by taking the derivative 
of the total energy. 



FIG. 13. (Color online) Kinetic energy of the particles vs. 
interaction strength for /3 = 0.2 and j3 = 10.2 for A = 50 par¬ 
ticles. The inset shows the interaction potential contribution 
to the total energy. The initial increase in interaction energy 
follows the increase in the interaction strength. However, be¬ 
yond a certain strength, the tendency of the fermions to avoid 
the impurity is more dominant. These plots are obtained by 
taking the derivative of the total energy with respect to c. Al¬ 
ternatively, the interaction potential energy is also obtained 
by using the two-particle correlation function at zero separa¬ 
tion. Both results are plotted in the inset showing remarkable 
agreement. 


V. CONCLUSION 

The problem of a single impurity interacting with a 
fermion background has attracted attention of the con¬ 
densed matter society for years. In this paper, we argue 
that an artificial gauge field coupling exclusively to the 
impurity is an effective tool to probe the physics of this 
system at any interaction strength. We consider a Eermi 
gas in a narrow ring trap and an artificial magnetic field 
coupling only to a single impurity. We solve this system 
exactly by using the BA for contact interactions and cal¬ 
culate the dependence of measurable quantities on the 
external magnetic flux. We observe this dependence for 
total energy, angular momentum of the charged particle, 
effective mass and the two-particle correlation function. 

Using an artificial magnetic field in this system has two 
advantages. The usual measurement tools such as expan¬ 
sion imaging become probes of thermodynamic quantities 
by comparing measurements at different flux values. For 
example, the change of the momentum carried by the im¬ 
purity caused by the magnetic field is a direct probe of the 
effective mass of the impurity. The second advantage is 
obtained by adiabatically increasing the flux value. Al¬ 
though the Hamiltonian of the system is periodic with 
flux, adiabatic evolution connects the ground state at 
zero flux to excited states at integer flux. In a cold atom 
experiment, such excited states can be expected to have 
long lifetimes due to total angular momentum conserva- 
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tion. Thus, we have calculated the physical properties 
for not only the ground state but also for excited states 
adiabatically connected to it. 

Our results show that the system can be described by 
a simple physical picture. The charged particle interact¬ 
ing with the background particles drags them along with 
itself around the ring. In the non-interacting limit, all of 
the angular momentum in the system is carried by the im¬ 
purity. As interactions are turned on, fermions which are 
close to the impurity in momentum are disturbed more 
and start to gain momentum. At the limit of infinitely 
repulsive interactions, the charged particle is effectively 
indistinguishable from the background and the total mo¬ 
mentum is shared equally between all particles. The pair 
correlation function also confirms this picture. The value 
of the correlation function near zero is mostly insensitive 
to the external flux while away from the correlation hole 
frequency of the Friedel oscillations sensitively depends 
on it. For strongly repulsive interactions, the effective 
mass saturates the total mass of the particles since it is 
dragging all of the background fermions along with itself 
around the ring. 

For attractive interactions, the impurity forms a bound 
pair with one of the fermions. The effect of dimer for¬ 


mation can be clearly seen in the angular momentum 
and the effective mass. For infinitely strong attractive 
interactions, angular momentum carried by the impurity 
saturates half the value of total angular momentum and 
the effective mass saturates twice the mass of the particle 
which confirm the presence of the dimer as a composite 
particle. 

The physical properties calculated in this paper are 
experimentally accessible through the standard tools of 
ultracold atom experiments. While artificial magnetic 
fields have been demonstrated in a variety of settings, 
they have not been used in combination with a toroidal 
trap to our knowledge. We believe our exact results 
would be relevant for such an experiment. 
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